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(Dated: March 10, 2008) 

We report the results of a search for a narrow resonance decaying into two photons in 1.1 fb _1 of 
data collected by the DO experiment at the Fermilab Tevatron Collider during the period 2002-2006. 
We find no evidence for such a resonance and set a lower limit on the mass of a fermiophobic Higgs 
boson of mhf > 100 GeV at the 95% C.L. This exclusion limit exceeds those obtained in previous 
searches at the Tevatron and covers a significant region of the parameter space B(hf — ► 77) vs. 
mh f which was not accessible at the CERN LEP Collider. 
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PACS numbers: 14.80.Ly, 12.60. Jv, 13.85.Rm 

In the standard model (SM), the Higgs field is respon- 
sible for both electroweak symmetry breaking and gen- 
crating elementary fermion masses. And while the SM 
describes our world at current experimentally accessible 
energies, the exact mechanism for electroweak symmetry 
breaking remains a mystery. 

Di-photon decays of the Higgs boson are suppressed 
at tree level, and in the SM such decays have a very 
small branching fraction, 10~ 3 — 1CP 4 . However, in a 
more general framework where the parameter content of 
the theory is richer, such decays can be enhanced. In the 
situation where the Higgs-fermion couplings are substan- 
tially suppressed, the full decay width of the Higgs bo- 
son would be shared mostly between the WW, ZZ, and 
77 decay modes. Such a scenario, the so-called "fermio- 
phobic" Higgs boson, arise in a variety of models, e.g. 
[TJ [2j [3] . In all these cases, for masses < 100 GeV, 
the Higgs boson dominantly decays to photon pairs. 

Experimental searches for fcrmiophobic Higgs bosons 
(h f ) at the CERN LEP Collider and the Fermilab Teva- 
tron Collider have yielded negative results. Mass lim- 
its have been set in a benchmark model that assumes 
that the coupling htVV (V = W ± ,Z) has the same 
strength as in the SM and that all fermion branching 
ratios (B) are exactly zero. Combination of results ob- 
tained by the LEP collaborations @J El [7] using the 
process e + e~ — > hfZ, hf — > 77 yielded the lower bound 
m h > 109.7 GeV at the 95% C.L. [8 . In Run I of the 
Tevatron, lower limits on mf lf from the DO and CDF col- 
laborations are respectively 78.5 GeV [S] and 82 GeV 
[TO] , using the processes qq' — > V* — > hfV, hf — > 77, 
with the dominant contribution coming from V = W . 

In this Letter we perform a search for the inclusive pro- 
duction of di-photon final states via the Higgsstrahlung 
and vector boson fusion processes: pp — > hfV — > 77 + X 
and pp — > VV — ► hf — ► 77 + X, respectively. The total 
integrated luminosity of the data used for this search is 
1.10 ± 0.07 fb _1 . 

The DO detector comprises a central tracking system in 
a 2 T superconducting solenoid, a liquid-argon/uranium 
sampling calorimeter, and a muon spectrometer. The 
calorimeter consists of a central section (CC) covering 
the pseudorapidity range \rj\ < 1.1, which is defined as 
r\ = — log[tan(|)] where 9 is the polar angle with respect 
to the proton beam direction, and two endcaps (EC) ex- 
tending coverage to |»7| < 4.2, each housed in a sepa- 
rate cryostat. The electromagnetic (EM) section of the 
calorimeter has four layers with longitudinal depths of 
2Xq, 2Xq, 7Xq, and 10Ao that provide full containment 
of EM particles (photons and electrons) . The calorimeter 
layers have transverse segmentation oi5(f>x5r] = 0.1 x 0.1 
(where <j) is the azimuthal angle), except in the third 
layer, where it is 0.05 x 0.05, which allows for accurate de- 



termination of the position of EM particles. Immediately 
before the inner layer of the central EM calorimeter there 
is a central preshower detector (CPS) formed of 2X$ of 
absorber followed by several layers of scintillating strips 
with embedded wavelength- shifting fibers. A complete 
description of the DO detector can be found in [TT] . 

We select events that satisfy single EM triggers which 
become fully efficient for EM showers with transverse 
momentum pt > 30 GeV. Photons and electrons are 
identified in two steps: the selection of EM clusters, and 
their subsequent separation into those caused by photons 
and those caused by electrons. EM clusters are selected 
from calorimeter clusters by requiring that (i) at least 
97% of the energy be deposited in the EM section of the 
calorimeter, (ii) the calorimeter isolation be less than 0.07 
(isolation is defined as [E tot (0A) - £ E m(0.2)]/£ , E m(0.2), 
where £'tot(0.4) is the total shower energy in a cone of 
radius R = ^(Ar/) 2 + (A0) 2 = 0.4, and E EM (0.2) is the 
EM energy in a cone with R = 0.2), (hi) the transverse, 
energy-weighted shower width be less than 0.04 rad (i.e. 
consistent with an EM shower profile), and (iv) the scalar 
Pt sum of all tracks originating from the primary ver- 
tex in an annulus of 0.05 < R < 0.4 around the cluster 
be less than 2 GeV. The cluster is then defined as an 
electron if there is a reconstructed track (or electron- like 
pattern of hits in the tracker) associated with it and a 
photon otherwise. We also consider EM jets (jets with 
a leading ir° or 77) defined as EM clusters that pass all 
cuts required for photon candidates except the track iso- 
lation requirement. We will refer to them as "j" or "jet". 
We select events that have at least two photons in the 
central calorimeter (\rj\ < 1.1) with transverse momenta 
Pt > 25 GeV. Events are required to have the primary 
vertex close to the beam axis and within 60 cm of the 
geometrical center of the detector. Identification of the 
primary vertex in the event is important, as it affects the 
calculation of the pt of a photon candidate and its track 
isolation. Despite the fact that photons do not leave 
tracks, the probability to reconstruct a primary vertex is 
high, 99.5%, due to the underlying event activity. 

The Higgs boson produced in the models considered 
has higher transverse momentum q'P than that of the 
two-photon system of the backgrounds. Therefore, we se- 
lect events with q T 7 > 35 GeV. For simplicity, we choose 
a fixed cut value which is below the optimal cut value for 
Higgs boson masses starting from 70 GeV. After all selec- 
tion criteria, we are left with 196 (1509) di-photon events 
with qip > 35 (qj? < 35) GeV for invariant masses above 
65 GeV. 

The dominant background comes from direct di- 
photon production (DDP) processes. The other major 
background comes from events in which jets are misiden- 
tified as photons: ~fj processes where a quark or a gluon 
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fragmented into an energetic 7r° or 77 and is reconstructed 
as a photon, and the multijet background where two jets 
are mis-identified as photons. 

Another source of di-photon background comes from 
events in which electrons are misidcntified as photons: 
the decay of a Z boson where electrons are reconstructed 
as photons if there are no associated tracks, and processes 
with one real electron coming from the decay of a W 
boson produced in association with a real photon or a 
jet misreconstructed as a photon. The veto of electron- 
like patterns of hits in the tracker reduces electron back- 
grounds by a factor of five, while keeping the photon 
efficiency high. We measure that (91 ± 3)% of photon 
candidates in Z/j* — > e + e~j data satisfy the anti-track 
activity requirement. The contribution of events with one 
or two real electrons is obtained by applying the prob- 
ability for an electron to fail the track requirement and 
be reconstructed as a photon (1.5i 1 ; s %) to the Z boson, 



This background 



-1.5' 

Drell-Yan, and W ± + X event yields, 
is estimated to be less than one event. 

We estimate the relative contributions of the 77, -fj, 
and jj backgrounds using the difference in the energy 
weighted width of the energy deposition in the CPS, 
&E ■ The width is generally narrower for photons than 
for jets. We construct one-dimensional templates as a 
function of x = <J^ PS for photons [G(x)] and jets 
The G(x) is constructed using radiative Z/j* — > 
(£ = e, /i) decays in data and the J(x) is taken from the jj 
data sample. From these we construct two-dimensional 
profiles for the three components 77, jj, and jj, as fol- 
lows: GG(x,y) = G{x) ■ G(y), GJ(x,y) = 0.5 • [G(x) ■ 
J(y) + J(x) ■ G(y)\, and JJ(x, y) = J(x) ■ J(y). Further, 
using these two-dimensional templates we construct a fit- 
ting function: cq ■ [GG(x, y) + c\- JJ{x, y) + C2 • GJ(x, y)]. 
The parameters are chosen so that cq is equal to the 
number of 77 events and responsible for the overall nor- 
malization, and Ci and C2 determine the contributions of 
jj and 7 j events relative to 77. 

For the di-photon candidate data sample, we make a 
two-dimensional distribution of <7§ PS . For each event we 
randomly decide whether the leading photon is plotted 
along the x- or the y-axis. We fit this distribution with 
the function defined above to determine the individual 
components: cq = 131 ± 22 ± 7 events, c\ = 0.35 ±0. 19 ± 
0.06, and C2 = 0.13 ± 0.28 ± 0.13, where the first error 
is the statistical error of the fit, and the second is the 
systematic uncertainty obtained from variations of the 
fitting range, binning of the templates, and the source of 
the photon template. 

The next step is to use the derived fractions to model 
the mass distribution of the di-photon candidate data. 
For this we need three mass templates: T 77 , T 7 j, and 
Tjj . We take T 77 from pythia MC [12] corrected for de- 
tector effects and reweighted with the K-factor derived 
from ResBos [T3] to account for the (next-to-) next-to- 
leading order, NLO (NNLO), effects. The other two tem- 



plates are taken from jj and jj samples, where we re- 
lax the calorimeter isolation, EM fraction, and energy- 
weighted shower width requirements in the definition 
of a jet in order to increase statistics in these tem- 
plates. We verify that relaxing the requirements do not 
alter the kinematics of the sample. We also correct the 
7j mass template for the admixture of jj events. We 
construct the background mass spectrum assuming the 
functional form N~ n ■ (T 77 + ci ■ Tjj + C2 • T~A where 
T!,,- are mass distributions normalized to 
one (see Fig. [I]), c\ and C2 are taken from the CPS fit 
above, and iV 77 is the expected number of DDP events 
from the MC. For the measured luminosity, we estimate 
iV 77 = 113±3.5(stat)±24(syst) events, which is in agree- 
ment with the c = 131±22±7 events derived from 
data. While these numbers, 113 and 131, are within the 
theoretical and experimental uncertainties, we choose to 
normalize the number of background events to the to- 
tal number of events observed in the data (normalization 
events are counted outside of the signal region, defined 
as a ±5 GeV window in diphoton mass centered at each 
hypothesized rrih s value) . By doing so we eliminate most 
of the background uncertainties, e.g. luminosity, renor- 
malization scale. 
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FIG. 1: Normalized distributions of the invariant mass, m, of 
77 (circles), jj (squares), and jj (triangles). 

Figure [2] shows the mass distributions in data with 
overlaid background predictions. The shaded regions cor- 
respond to the expected background error bands. The 
inner band represents the statistical uncertainty of the 
mass templates, while the outer corresponds to the sys- 
tematics due to variation in the one-dimensional cr§ PS 
templates. We assign an additional 100% uncertainty 
that includes any possible change in the shape of the 
mass templates due to the relaxed definition of a jet. 

Signal events are generated for a range of mass points 
from 70 GeV to 150 GeV in 10 GeV steps. We use the 
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FIG. 2: Di-photon mass distribution of the data (squares) 
with the overlaid background prediction (triangles), and the 
expected signal distribution for mhf = 100 GeV in the bench- 
mark model. 



pythia event generator followed by a detailed GEANT- 
based [H] simulation of the DO detector. The signal ef- 
ficiencies, e slgnal , are derived from the MC. Table [i] lists 
signal efficiencies after correction for trigger inefficiency 
and scaling by the ratio of efficiencies in data and MC 
(« 95% per photon) obtained from the electron recon- 
struction efficiency in Z — » e + e~ events. Note that the 
photon requirements are chosen in such a way that the 
MC correctly reproduces differences between electrons 
and photons as confirmed in Z — > e + e _ 7 events. Ta- 
ble |T] also shows the number of observed di-photon can- 
didate events in data in 10 GeV mass windows and the 
corresponding background estimates with associated un- 
certainties. The width of the mass peak is dominated by 
the detector resolution and varies between 2.8 GeV and 
5.2 GeV. The size of the optimal mass window varies be- 
tween 8 GeV and 15 GeV, but for simplicity we use a 
fixed value of 10 GeV. The acceptance of the mass win- 
dow cuts varies between 94% and 66% for mh f — 70 — 150 
GeV. In the same table we provide the theoretical bench- 
mark branching ratio, B(h — > 77) (15j . and the NLO 
cross section, er^ LO , for the sum of the signal processes 
pp -> VV 

v2hv na. 



hf and pp — ► hfV obtained with vv2h and 



We perform a counting experiment in the 10 GeV mass 
windows, and in the absence of an excess of di-photon 
events, we set an upper limit on the product of the Higgs 
boson production cross section and di-photon branching 
ratio <Jh f ■ B(hf — ► 77) at 95% C.L. Limits are calculated 
using the modified frequentist CLg method [17]. Table [I] 
shows the expected and observed limits. The present 
study excludes fermiophobic Higgs bosons of mass up to 
100 GeV at the 95% C.L. This is the most stringent limit 



to date at a hadron collider. In Fig. [3] we present our 
results as limits on the branching ratio in the parame- 
ter space B(hf — > 77) vs. m/,, obtained from a ratio 
of the above limits and cr^ L °. The regions above the 
experimental points correspond to the excluded values 
of the branching ratio. This study significantly improves 
the LEP limits at intermediate mass values, e.g. by more 
than a factor of four at mh, — 120 GeV, and extends sen- 
sitivity into the region not accessible at LEP, m/,, > 130 
GeV. 
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FIG. 3: B(hf — > 77) limits as a function of the Higgs mass. 
The theoretical BQif — > 77) curve for the benchmark model 
as well as the observed B(hf — > 77) limits from DO Run I 
and LEP are overlaid. The shaded regions correspond to the 
excluded values of the branching ratio. 
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a(pp — > hf + X) ■ BQlf — ► 77) (pb) 
mk, (GeV) data background e slgnal (%) expected limit observed limit Run I limit cr^ LO (pb) B(hf —* 77) 



70 


35 


24.5 ± 4.6 


6.9 ± 0.5 


0.15 


0.29 


0.46 


1.5 


0.81 


SO 


33 


27.2 ±5.0 


7.9 ±0.6 


0.14 


0.20 


0.44 


1.0 


0.70 


90 


24 


27.4 ± 5.4 


9.8 ±0.8 


0.11 


0.089 


0.37 


0.75 


0.41 


100 


24 


23.7 ±4.8 


10.3 ±0.8 


0.10 


0.10 


0.35 


0.55 


0.18 


110 


14 


17.7 ±4.4 


11.2 ±0.9 


0.085 


0.061 


0.34 


0.42 


0.062 


120 


11 


13.4 ±3.7 


11.3 ±0.9 


0.070 


0.058 


0.33 


0.32 


0.028 


130 


9 


11.7 ±3.3 


11.2 ±0.9 


0.065 


0.053 


0.33 


0.25 


0.019 


140 


8 


9.5 ±2.8 


11.7 ±0.9 


0.058 


0.052 


0.32 


0.19 


0.0061 


150 


12 


6.3 ± 2.1 


11.7 ±0.9 


0.051 


0.10 


0.32 


0.15 


0.0020 



TABLE I: Input data for limit calculation and 95% C.L. limits on cross section times branching fraction. Quoted are the total 
uncertainties that are used in the limit calculation. 
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